THE SOLUTION OF THE A r T-SYSTEM FOR ARBITRARY 

BOUNDARY 



PHILIPPE DI FRANCESCO 

Abstract. We present an explicit solution of the A r T-system for arbitrary boundary 
conditions. For each boundary, this is done by constructing a network, i.e. a graph with 
positively weighted edges, and the solution is expressed as the partition function for a 
family of non-intersecting paths on the network. This proves in particular the positive 
Laurent property, namely that the solutions are all Laurent polynomials of the initial data 
with non-negative integer coefficients. 
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I. Introduction 

In this paper we study the solutions of the A r T-system, namely the following coupled 
system of recursion relations for a,j, fceZ: 

(1-1) T a jjfc+lTci ,j,k— 1 r ^a,j+\,k^a,j— l,k Ti+lj^fcTi— l,j,k 

for a £ I r = {1,2, ...,r}, and subject to the boundary conditions 
(1-2) T j,k = T r +ij,k = 1 (j, k £ Z) 

This system arose in many different contexts. The system (jl.ip and its generalizations 
were introduced as the set of relations satisfied by the eigenvalues of the fused transfer 
matrices of generalized quantum spin chains based on any simply- laced Lie algebra g [2] [15] ; 
in this paper we restrict ourselves to the case q = sl r+ %, but we believe our constructions 
can be adapted to other g's as well. 

With the additional condition that T Qj0 ,fc = 1, A; £ Z and the restriction to j £ Z + , the 
solutions of fll.mi.2p were also interpreted as the q-characters of some representations of 
the affine Lie algebra U q (sl r+ i), the so-called Kirillov-Reshetikhin modules [UJ, indexed by 
a £ I r = {1, 2, r} and j £ Z + , while k stands for a discrete spectral parameter [TT] . 

The same equations appeared in the context of enumeration of domino tilings of plane 
domains J2U], and was studied in its own right under the name of octahedron equation [13] 
[13J. As noted by many authors, this equation may also be viewed as a particular case 
of Pliicker relations when all T's are expressed as determinants involving only the Ti^'s. 
These particular Pliicker relations are also known as the Desnanot-Jacobi relation, used by 
Dodgson to devise his famous algorithm for the computation of determinants [8]. In [19] . 
this equation was slightly deformed by introducing a parameter A before the second term 
on the r.h.s. and used to define the "lambda-determinant", with a remarkable expansion on 
alternating sign matrices, generalizing the usual determinant expansion over permutations. 
Here we will not consider such a deformation, although we believe our constructions can be 
adapted to include this well (see J2U] for a general discussion, which however does 

not cover the A r case). 

Viewing the system (11. ip as a three-term recursion relation in k £ Z, it is clear that the 
solution is entirely determined in terms of some initial data that covers two consecutive 
values of k, say k = 0, 1 and say all j £ Z. In [6], an explicit expression for T a j^ was 
derived as a function of the initial data x = {T a j i0 , T at j t i} ae i rj j e z- It involved expressing 
first Tij k as the partition function for weighted paths on some particular target graph, with 
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weights that are monomials of the initial data, and then interpreting T a j t k as the partition 
of non-intersecting families of such paths. This interpretation was then extended to other 
initial data of the form 

(1.3) x k = {T a j :ka , T a j tka+ i} aeIr j e z 

where k = (k±, ki, • k r ) G Z r is a Motzkin path of length r — 1, namely k a+ \ — k a G 
{0,1,-1} for all a = l,2,...,r — 1. In this construction, for each Motzkin path k, the 
expressions for the T a j^ in terms of the initial data x^ are also partition functions of 
weighted paths on some target graph T^. 

The equation (11.11) is also connected to cluster algebras. In Ref. [3], it was shown that 
the initial data sets x k form a particular subset of clusters in a suitably defined cluster 
algebra. Roughly speaking, a cluster algebra [10] is a dynamical system expressing the 
evolution of some initial data set (cluster), with the built-in property that any evolved 
data is expressible as Laurent polynomials of any other data set. This Laurent property 
or Laurent phenomenon turned out to be even more powerful than expected, as all the 
known examples show that these polynomials have non-negative integer coefficients. The 
positivity conjecture of [TU] states that this property holds in general. As an example, the 
above-mentioned lambda-determinant relation may be viewed as an evolution equation in 
the same cluster algebra as in [3], with A as a coefficient: the existence of an expansion 
formula of the lambda-determinant on alternating sign matrices is a manifestation of the 
positive Laurent phenomenon. As another example, the explicit expressions of [6] for the 
solutions T a j t k of the A r T-system as partition functions for positively weighted paths gives 
a direct proof of Laurent positivity for the relevant clusters. 

However, the set of initial data x^ (II. 3p covered in [B] is limited to sets of T a j^s with 
fixed values of k = k a independently of j. The most general set of initial data should also 
allow for inhomogeneities in j, namely values of k = k a> j varying with j as well. It is easy 
to see that the most general boundary condition consists in assigning fixed positive values 
(daj)aei r ;jez to T a> j } k aJ along a "stepped surface" (also called solid-on-solid interface in the 
physics literature), namely such that \k a+ ij — k a j \ = 1 and \k a j + i —k a j\ = 1 for all a G I r 
and j G Z. 

In this paper we address the most general case of initial data for the A r T-system (11.11 
II. 2p . As we will show, initial data are in bijection with configurations of the six- vertex 
model with face labels on a strip of square lattice of height r — 1 and infinite width. For 
any such given set of initial data, we derive an explicit expression for the solution T a 
(ll.mi.2p as the partition function for a non-intersecting paths on a suitable network, in 
the spirit of Refs. |E] and [18] , and with step weights that are Laurent monomials of the 
initial data. This completes the proof of the Laurent positivity of the solutions of the A r 
T-system for arbitrary initial data. 

The paper is organized as follows. 

Our construction was originally inspired by Ref.[T] which basically deals with the case 
of A\ under the name of "frises": the latter is reviewed in Section 2, where we make in 
particular the connection between the "frise" language and the solutions of the A\ T-system 
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with arbitrary boundary data. Roughly speaking, the solution is expressed as the element 
of a matrix product taken along the boundary. 

A warmup generalization to the case of A 2 is presented in Section 3, with the main 
Theorem 13.41 giving an explicit solution for arbitrary boundary data, also as an element of 
a matrix product taken along the boundary. 

Section 4 is devoted to the general A r case. Starting from the path solution of [6] for some 
particular initial data, we construct various transfer matrices associated to the boundary, 
with simple transformations under local elementary changes of the boundary (mutations). 
For convenience, boundaries are expressed as configurations of the six-vertex model in an 
infinite strip of finite height r — 1. These in turn encode a network, entirely determined by 
the boundary data. The final result is an explicit formula Theorem 14. 121 for the solution of 
the A r T-system as the partition function for families of non-intersecting paths. 

In Section 5, we study the restrictions of our results to the Q-system. 

A few concluding remarks are gathered in Section 6. 

2. At T-SYSTEM AND FRISES 

In this section, we first review the results of |lj, and then rephrase them in terms of 
solutions to the A% T-system for arbitrary boundary conditions. 

2.1. Frises. 

2.1.1. Frist equation. The frise equation readfl 

(2.1) Ma+l,6-lM a ,6 = 1 + Wo+l,6«o,6-l 

for a, b G Z. 

2.1.2. Boundaries. The most general (infinite) boundary condition is along a "staircase", 
made of horizontal (h) and vertical (v) steps of the form h : (x,y) — > {x + l,y) and 
v : (x,y) —> (x, y + 1), giving rise to a sequence of vertices (xj,yj), j G Z. To each vertex 
of the sequence we attach a positive number aj, j G Z, and the boundary condition for the 
system ( 12. ip reads: 

(2.2) u Xj>yj = aj (j G Z) 

The simplest such boundary is the sequence ...hvhvhv..., say with variable a 2x at vertex 
(x, x) and c^+i at vertex (x, x + 1), x G Z. We refer to it as the basic staircase boundary. 
The problem is now to find the solution u a ^b of (12. ip with the boundary condition ( 12. 2p . 

2.1.3. Projection of (x,y) on the boundary and step matrices. The general solution at a 
point (x, y) to the right of the boundary is expressed solely in terms of the values (I2.2p 
taken by u along the "projection" of (x, y) onto the boundary, defined as follows. 



Our convention corresponds to b 



—b in those of Ref pQ. 
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FIGURE 1. A typical boundary for the frise and the projection of a point (x, y) 
onto it. The vertices of the projection on the boundary are represented with blue 
circles. Here, the word w(x,y) reads v 2 hvhvh 3 . 

Definition 2.1 (Projection). The projection of (x,y) onto the boundary {(xj,yj)}j e z is 
the sequence (xj,yj), j = t, t + 1, ...,t', where y t = y, x t > = x, and the first step t — > t + 1 
is vertical, while the last step t' — 1 — > t' is horizontal. 

This is illustrated in FigJTJ Alternatively the projection of (x, y) is coded by the word 
w(x,y) =v...h of length if — t with letters h and v, starting with v and ending with h 
and coding the succesion of horizontal (h) and vertical (v) steps along the boundary be- 
tween (x t ,yt) and (x t >,yt>)- We also define the corresponding sequence of boundary weights 
a.(x,y) = (a t ,a t+ i, ...,a#). 

Definition 2.2 (Step matrices). We define the two horizontal and vertical matrices 

(2.3) *M) = £(1 °) "M) = i(S i) 

2.1.4. Solution. Given some boundary conditions, we associate to the word w and the 
sequence a the following 2x2 matrix product 

(2.4) M(w, a) = V(a t , a t+1 ) ■ ■ ■ H(a t ,_ u a t >) 

where the product extends over all the intermediate steps % — > i + 1 between t and t' as 
coded by w, and involves the matrix H(ai, a i+1 ) if the step i — > i + 1 is h and V(ai, a i+ i) if 
it is v. The result of [1] takes the following form: 

Theorem 2.3 ([I]). The solution of (12. ip subject to the boundary condition (12. 2p reads: 

(2.5) u XiV = a t > (M(w(x, y), a(x, y))) hl 
with M as in ([23]) . 
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All matrices V, H having elements that are positive Laurent monomials of the initial 
data, the general Laurent positivity of the solution follows: 

Corollary 2.4. The general solution of ( 12. ip subject to the boundary condition (12. 2p is a 
Laurent polynomial of its initial data {aj}j(zz, with non-negative integer coefficients. 

Example 2.5 (The basic staircase boundary). For any x > y G Z> , we have a projection 
on the boundary with w(x,y) = {vh) x ~ y , and a.(x,y) = (a 2y , a 2y +i, ...,a 2x ). We deduce that 

x-l 

(2.6) M(w(x, y), &(x, y)) = JJ V(a 2i , a 2i+1 )H(a 2i+1 , a 2i+2 ) 

i=y 

and the solution reads: 

(2.7) u X)V = a 2x (M(w(x, y), &(x, y))) hl 
Explicitly, we compute the two-step matrix: 

(2.8) M(a, b, c) = V(a, b)H(b, c) = l - {^f j) 

2.1.5. Mutations. Note that we may move from one boundary to another by elementary 
"mutations"[l, namely the local substitution (v , h) — > (h, v) on the boundary (forward 
mutation) or (h, v) — > (v,h) (backward mutation), while the sequence a is updated using 
the frise relation (12. ip . In particular, we may in principle reach any boundary from the 
basic staircase one, by possibly infinitely many such mutations. 

The effect of such a mutation is easily obtained by computing the corresponding matrix 
transformation within M(u>,a). It basically corresponds to the following identity: 

Lemma 2.6. For all a,b,c> 0, we have: 
b - c 

(2.9) \ V{a, b)H{b, c) = H(a,x)V(x,c), x = 




This may be understood as a matrix representation of the mutation via the commutation 
of the matrices V and H, which acquire the new boundary value x in replacement for b. 
This mutation affects all values of u m ^ p such that the projection of (m,p) contains the new 
boundary point with value x. 

We may deduce the general formula (I2.5P from that for the basic staircase boundary, by 
induction under mutation. In general a mutation simply switches two consecutive matrices 
VH — > HV in the product M. We must be careful with mutations that update the extremal 
vertices of the projection of (x,y), namely in the two cases: (i) when w starts with vh, 
updated into hv or (ii) when w ends up with vh, updated into hv. We note however that 
H(a,b)ij = bV(a,b)j t i = S^j hence in the updated matrix M we may: (i) drop the first 



2 The term "mutation" is borrowed from cluster algebras, as this elementary move indeed corresponds 
to a mutation in the associated cluster algebra of [5] . 
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matrix factor H (ii) drop the last matrix factor V, but replace the scalar prefactor by the 
new updated vertex value, and the formula ( 12 .5p follows. 

2.2. The Ay T-system. 

2.2.1. T-system. The Ay T-system reads: 

(2.10) Tj^+yTj^k-i — Tj + y.kTj-y t k + 1 

where we use the shorthand notation Tj k = Tyj jk for j, k G Z. Note that this splits into 
two independent systems for fixed value of j + k modulo 2. 

In the case when j + k = modulo 2, we immediately see that changing to "light cone" 
coordinates: a = and b = we have that u a> b = Tj >k _y satisfies the frise equation 
12. ip . So the two problems are equivalent. Analogously, when j + k = 1 modulo 2, we take 



a = 1 and b = J and u a ^ = Tj >k . 

2.2.2. Boundaries. The fundamental boundary for the T-system is obtained by fixing the 
values of say T 2 j + y t o and T 2 j t y for all j 6 Z. It corresponds to the basic staircase boundary 
in the case j + k = 1 modulo 2 above, with T 2 ^y = a 2 j and T 2 j +lfi = a 2 j +1 . 

Other boundaries are mapped in an obvious manner. 

2.2.3. Path solution. In [6], an explicit path formulation was derived for the solution Tj t k 
for the fundamental boundary condition. Defining the 4x4 transfer matrix 



7(u, v, w) 



(0 


1 





o\ 


u 





1 








V 





1 







w 


V 



we have 



Theorem 2.7 (|6J). The solution of the Ay T-system (I2.10p for the fundamental boundary 
condition with initial data {T 2 j+y t o, T 2 j > y}jez reads for j + k = 1 mod 2: 

/j+k-y 

T j>k = T j+k>0 I J| 7(ui,Vi,Wi) 

\i=j-k / l l 

m 1 1\ Ti,y 1 1 

I 2 - 11 ) Ui = 7f > Vi = — — , Wi = — 

J-i+1,0 J-iflJ-i+l,! Ui 

Here the matrix 7(i, % + 1) = 7(ui, Vi,Wi) is interpreted as the transfer matrix from time 
i to time i + 1 for weighted paths with steps a — » a ± 1 on the integer segment [0, 3], with 
time- dependent step weights — > 1 : u^, 1 — > 2 : and 2 — > 3 : Wi, the other weights being 
equal to 1. 
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2.2.4. Gauge invariance. The above formula (12.111) remains clearly unchanged if we trans- 
form the matrix T into the following: 7(i, i + 1) — > 7{i, i + 1) = Li7(i, i + 1)L^ 1 , for any 
invertible matrix Lj such that (Li)xj = and = 6ji. 

To make the contact with the frise solution, let us define T as above, by use of the matrix 

(\ \ 

10 

T tfi 

\0 T iA J 

2.2.5. Comparison with the frise solution. Let us now compute the "two-step" transfer 
matrix at times i, i + 1 for i = j + k modulo 2: 7(i, i + 2) = T(z, z + l)7(i + 1, i + 2), with 
the weights as in (12. lip , namely: 



U: 



h+x 



V i+ X 



Wi = — 
Ui 



1 



1 



w 



i+1 



a i+2 °i+xbi+2 

where we have introduced a { = T ij0 , a i+1 = T i+l ,x, ai+2 = T i+2> o, h = T itl , b i+1 = T i+lfi , 
and b i+2 = T i+2 ,x, while U = diag(l, 1, a h bt), L i+1 = diag(l, 1, b i+l , a i+1 ) and L i+2 = 
diag(l, l,a i+ 2,bi +2 ). 



We get: 



7{i,i + 2) 



a.i+2 

1 

a.i+2 

V 



o 

X+bjbj+2 
bi+ibi+2 


1 

h+2 



1 


\ 





1 


bi+2 


X+aiO.i+2 





ai+ia.i+2 





bi+2 / 



This 4x4 matrix clearly decomposes into two independent 2x2 linear operators acting 
respectively on components 1,3 and 2,4. The corresponding matrices are respectively: 

i _ \ 

H(ai, a i+ x)V(a i+ x } a i+2 ) 



T(M + 2) 



a.i+2 
X+ajai+2 

^ai+2 ai+idi+2 
'X+bibi+2 1 




V{hM^)H{b i+ x,b 



i+2) 



We may now use the gauge-transformed and reduced two-step transfer matrix 7(i, i + 2) 
instead of 7 in (12. lip . Indeed, in the product over steps from j — k to j + k — 1, we may pair 
up consecutive 7 matrices in (12. lip to express it in terms of the CP's, and then substitute 
the latter with the CP's, leading to: 

/fc-i 

T hk = T j+kfi Yl Hj -k + 2i,j-k + 2i + 2) 



1,1 
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Noting moreover that if (a, b)x,j = bV(a, b)j t i = 8j t i, we may rewrite this as 

(k-2 
n v ^ 
-fe+2i+l) a j-k+2i+2)H(aj^ k+2 i + 2, Oj^k+2i+3) 
i=0 / 1,1 

Assuming that j + k — 1 modulo 2, we see that in light-cone coordinates with x = 

and y = , equation (12 . 12 j) amounts to equations fl2.6H2.7p . as the projection of (x,y) 

on the boundary staircase starts at t = j — k + 1 and ends at t' — j + k — 1. 

2.2.6. Mutations and arbitrary boundary. As in the frise case, this identification gives us 
access to mutations, via the VH -h- HV identity (12. 9p . Starting from (I2.12p . we may 
iteratively apply forward/backward mutations to the basic staircase boundary to get any 
other boundary (up to global translations) of the form {Tj^Ajez with a sequence kj G Z 
such that |fcj+i — fcj| = 1. Let us denote by (ja,kj ) and (ji,^) the extremities of the 
projection of (j, fc) onto the boundary, namely such that j — kj — j — k, ji + kj ± = j + k, 
jo maximal and ji minimal. 

We deduce that the general solution for arbitrary staircase boundary reads: 

Tj,k = Tj lt k h {v{T j0jkjQl T jQ+ i^ jo+1 )...H(T h _i^ h _ 1 ,T jl ^ Ji ^ ^ 

where the product is taken along the projection of (j, k) on the boundary, with a matrix V 
per vertical step and H per horizontal step. 

3. The A 2 T-system with arbitrary boundary 
Before going to the general A r case, we derive the A 2 solution in detail. 

3.1. T-system. The A 2 T-system reads: 

Tij } k+iTij,k-i — Tij+i } kTij-i } k + T 2 j j k 
(3.1) T 2 j^+iT 2 j t k-i = T 2 j +1^2, j-i,k + T\j t k 

for j, fc 6 Z. Note that this splits again into two independent systems for T a j fc with fixed 
value of a + j + k modulo 2. These indices run over two consecutive layers of the centered 
cubic lattice a = 1 and a = 2, which form two square lattices, the vertices of the second 
layer lying at the vertical of the centers of the faces of the first layer. 

3.2. Boundaries. The fundamental boundary considered in Ref. [6] involves fixing the 
values of the T QJjfc with a = 1,2 k = 0,1, and j G Z, with fixed parity of a + j + k (say 
even). We refer to this boundary as the basic staircase boundary, in reference to the A\ case. 
It can be viewed as an infinite strip made of a succession of four kinds of vertices (see FigfS]). 
We may also view this strip as a succession of edges of the form ej = (1, j, 0) — (2, j, 1), 
f J+1 = (l,j + l,l)-(2,j + l,0),e, +2 = (l,j + 2,0)-(2,j + 2,l),etc. for j - 1 G 2Z (thick 
black lines in Figj2]). Two such consecutive edges define a tetrahedron. The basic staircase 
may therefore be viewed as the alternating succession of two kinds of tetrahedrons denoted 
by A (defined by ej, fj+i) and B (defined by fj-i, ej). 
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(2,j + 3,0) (2,j + 4,1) 




(1,1-1,1) 

FIGURE 2. The basic staircase boundary for the A 2 T-system. We have indi- 
cated the corresponding succession of edges (thick black line) and the two types 
of tetrahedrons A,B that connect them. 



3.3. Solution for the basic staircase boundary. In Ref. [6], the solution Tij t k was 
expressed in terms of paths on a target graph with 6 vertices and with time-dependent 
edge weights involving only the boundary values. These weights are coded by a 6 x 6 
transfer matrix. Defining: 



(3.2) 



7(s, t, u, v, w) 





1 













s 





1 














* 





1 


1 











u 

















V 








1 


Vo 











w 


0/ 



and using the notation 
(3.3) 

^1,U , ^2,i,l 



-^1,1+1,0^2,1-1, 



i+1,0 



22,1,0^2,1+1,1 



1 l.i+1,0 1 l,i,0 J l,i+l,l i l,i,l J 2,i,0 

we have: 

Theorem 3.1 Q6J). The solution of the A 2 T-system for a = 1 reads: 



(3.4) 



/j+k-i 

Ti,j,k = T ltj+ki0 I j j 7(s i ,t h u h v i ,w i ) 

\i=j-k / 1 



T 



2,i+l,0 
T2AA 



3.4. Reduced transfer matrix. As before we note that the two-step transfer matrix 
7(i, i+2) = 7(i, i+l)T(z+l, i+2), with 7(i, = 7(si,ti, Ui, V{, Wi), is again decomposable 
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into two linear operators acting respectively on components (1, 3, 6) and (2, 4, 5). Explicitly: 



T(s, t, u, v, w)7(s ,t ,u ,v ,w) 








1 













s + t' 





1 


1 





ts' 





t + u' + v' 








1 





ut' 





;/ 


u 








vt' 





V 


v + w' 





Vo 





wv' 








w) 



Defining: 

/ s l+1 1 

i + 2) = Us i+1 U + u i+1 + v l+ i 1 

y WiV i+ i Wi 

we may rewrite ( 13. 4p as: 

k-l 

(3.5) T 1Jtk = T ld+kfi ( Yl ?{j - k + 2i, j - k + 2i + 2)) 

i=0 



3.5. Gauge transformation, tetrahedron decomposition, and frise analogy. As 

before, we note that any gauge transformation of the form i + 2) — > i + 2) = 
LjT(i,i + 2)L~ f 1 2 and such that (Li)ij = = leaves the formula (13 .5p invariant. 

We choose 



Defining 



we have 



Li 




















1 







M(a, 6, c, m, v, w) 



u u I an I a _a_ 

c be fry t;u> «> 



o 1 



i + 2) — M (Tijfi, Ti^+i^, ^~2,i,l, ^2,2+1,0; ^2,i+2,l) 

The matrix M may be further decomposed as follows: 



M(a, 6, c, m, v, w) = A(a, b, u, v)B(b, c, v, w) 
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where 




A(a, b, u, v) 



1 0' 

u au a 

b bv v 

1 



'b i o 

c c 



B(b,c,v,w) =010 



,0 1 " 



W W I 



In view of our interpretation of the boundary (see Figj2]), the matrices A and B may be 
associated to the tetrahedrons A and B. The arguments are the values of T at the vertices 
of the tetrahedrons. We write 

A(i, i+l) = A(T 1>it i, T hi+lj0 , T 2AQ , T 2 , i+ i 7 i) and B(i, = B(T 1A0 , T 1>i+1>1 , T 2iiil , T 2>i+ i j0 ) 
and finally we may rewrite (I3.5P as: 

k-l 

(3.6) T ldtk = T 1J+kfi (j[A(j-k+ 2i, j - k + 2i + l)B(j - k + 2i + 1, j - k + 2i + 2)) 

i=0 

We may now interpret this result as a generalization of the frise result. By analogy with 
the frise solution, let us define the projection of (1, j, k) on the boundary as the portion of 
the boundary between the edge fj-k+i and the edge fj + k-i- We have: 

Theorem 3.2. The general solution of the A 2 T-system with the basic staircase boundary 
for a = 1 reads: 

fc-2 

(3.7) T ldtk = Tij +k ^i(jjB(j-k+2i+l,j-k+2i+2)A(j-k+2i+2,j-k+2i+3) S ) 

i=0 ' 

Proof. We start from (I3.6P and use the fact that A(a,b,u,v)ij = bB(a,b,u,v)j t i = 5j y i to 
eliminate the first (A) and last (B) matrices in the product on the r.h.s. □ 

The product extends over the sequence of tetrahedrons along the projection of (l,j,k) 
onto the boundary. We may think of the two tetrahedron matrices A, B as a generalization 
of the horizontal and vertical matrices of the A x case, but more general boundaries involve 
four more such matrices, as discussed below. 



3.6. Other boundaries: two tetrahedrons and four parallelograms. The most gen- 
eral boundaries are obtained from the basic staircase via local elementary moves (for- 
ward/backward mutations) corresponding to one application of one of the system relations. 
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(2,j+3,0) (2,1+4,1) 




(IJ-1,1) (1,1,2) (1,j+1,3) 



FIGURE 3. Various mutations of the basic staircase boundary for the A2 1 
system. We have indicated the tetrahedrons A,B and the parallelograms C,D,E, 
that may connect two consecutive edges of the boundary. 
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a b 



FIGURE 4. Triangle decomposition Ai,A2,Bi,B2 of the tetrahedrons A,B and 
parallelograms C,D,E,F. We use the same letter for different triangles that share 
the same transfer matrices (see eq. (|3.10j )). 



The effect is of flipping a single thick edge of the boundary in the following manner: de- 
noting by ej t k = (1, j, k) — (2,j, k + 1) and fj^ = (1, j, k + 1) — (2, j, k), we have the two 
possible elementary (forward) moves: 

(3.8) . . . , e^fc, fj t k+i ■ ■ ■ 

(3.9) . . . , fj t i~, ej jfc+ i . . . 

It is easy to see that this gives rise to six possible relative positions for two consecutive 
edges: 

( e i,fc; fj+i,k), [ej,k, ej+i^+i), (e^k, ej+i,k-i), (fj,k, e^+i^), (fj,k, fj+i,k+i), (fj,k, fj+i,k-i) 

which define two tetrahedrons and four parallelograms, respectively denoted by A,C,D,B,E,F 
(see Fig|3J. 

Note that the boundary is entirely specified by a Motzkin path and the edge at one of 
its vertices. Indeed the transition from an edge to the next changes k — > k, k + I or k — 1, 
and the nature of the edge (e or /) is switched only if k is unchanged. So keeping a record 
of the variable k is sufficient, and the record is a Motzkin path (or equivalently an infinite 
word in 3 letters). 



3.7. Triangle decomposition. As already mentioned, we may further decompose the 
tetrahedrons A and B, as well as the parallelograms C,D,E,F, into pairs of triangles, as 
indicated in FigJH To all triangles labelled Ai,A2,Bi,B2, we associate the following 3x3 
"triangle" matrices with 3 parameters equal to the values of T a j tk at their vertices. We 
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have: 

A o o\ 

A 1 (a,b,u) = H | o A 2 (&,u,«) 
V° V 

/! I o\ 

(3.10) Bi(a,6,«) =010 B 2 {b,u,v) 

\0 1/ 

The two tetrahedrons A,B, correspond to the matrices 

A(a,b,u,v) = Ai(a,b,u)A 2 (b,u,v) = A 2 (a,u,v)A%(a,b,v) 

(3.11) B(a,b,u,v) = Bi(a, b, u)B 2 (b, u, v ) — B 2 (a, u, v)Bi(a, b, v ) 

independent of the two triangle decompositions. The parallelograms C,D,E,F of FigJUhave 
unique triangle decompositions, to which we attach the following matrices: 

C(a,b,u,v) = Ai(a, b, u)B 2 {b^ u, v) - 
D(a,b,u,v) = A 2 (a,u,v)Bx(a,b,v) - 
E(a,b,u,v) = B 2 (a,u,v)Ai(a,b,v) - 

F(a,b,u,v) = Bi(a,b,u)A 2 (b,u,v) = 

Note that the products are taken in a specific order, namely that the matrix for the triangle 
which lies on the left of the diagonal of the parallelogram multiplies that on the right from 
the left. 




3.8. Mutations via triangles and the general formula. The general formula for Ti j t h 
for an arbitrary boundary reads as follows. First, we may view the boundary above in yet 
another manner, by projecting it vertically onto the bottom plane, as illustrated in Fig. 
It is the superimposition of the two sets of boundary vertices in the bottom (resp. top) 
layers (represented as filled (resp. empty) circles), which are both staircase boundaries of 
the type considered in the A\ case, drawn on the two corresponding shifted square lattice 
layers in thick (resp. dashed) lines. These two staircases are constrained by the condition 
that their vertices must be connected via e or / edges only (diagonal thick black lines in 
Fig E]). 
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, >. s 

: (Lie's) (i,j,k) 



(a) (b) 

FIGURE 5. The projection of a point (l,j, k) onto a typical boundary for the A 2 
T-system (a), viewed in vertical projection onto the bottom plane. The vertices of 
the bottom (rep. top) layers are represented as filled (resp. empty) circles. The 
boundary edges are represented as thick diagonal lines. We have indicated (b) the 
two types of forward mutations corresponding to e — > f and / — > e as in (J3.8H3.9p . 

Definition 3.3. The projection of the point (l,j,k) onto the boundary is the portion of 
boundary (finite sequence of edges) between the edges containing (l,j ,k ) and (l,ji,ki), 
repectively such that j — k = j — k and j\ + k\ = j + k with j Q maximal and j\ minimal. 

This coincides with the definition for the A\ T-system, using the staircase of the bot- 
tom layer only. The corresponding finite sequence of edges corresponds alternatively to 
a finite sequence of triangles according to the decomposition above, modulo the two-fold 
ambiguities of decomposition of the tetrahedrons A and B. 

To this sequence we associate the matrix: 

(3.12) M(j,k)= J] Z(x,y,z) 

triangle Z 
w/ vertex values x,y,z 

where for each triangle Z=A 1 ,A 2 ,B 1 ,B 2 , along the sequence we multiply by the correspond- 
ing triangle matrix Z = Ai, A 2 , Bi, B 2 . Note that, due to the identities ( 13. lip , this def- 
inition is independent of the particular choice of triangle decomposition of the possible 
tetrahedrons along the boundary. We have: 

Theorem 3.4. The solution of the A 2 T-system for arbitrary boundary reads for a = 1: 



(3.13) 



Ti,j, k = Ti JlM M(j, k) ltl 
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with M(j,k) as in ( 13. 121) with the product extending over the projection of (l,j,k) onto 
the boundary. 

Before proving the Theorem by induction under mutation, let us describe the mutations 
of the boundary in more detail. The two possible mutations fl3.8H3.9p correspond to a local 
transformation of the chain of triangles that forms the boundary, namely it replaces a pair of 
adjacent triangles sharing the initial boundary edge with a new pair of adjacent triangles 
sharing the mutated boundary edge. Using the definition (13.101) . we get the following 
Lemma, generalizing Lemma 12.61 

Lemma 3.5. For all a, b, c, u, v, w > we have: 
A, 

qq -|- 

Bi(d) 6, u)Ai(b J c, u) = A 1 (a J x, u)Bi{x, c,u), x = — - — 



uw ~\- b 

A 2 {b,u,v)B 2 {b,v,w) = B 2 (b,u,y)A 2 (b,y,w), y 



A V - 7 — 7 / — & \ 7 7 — / ■£> \ 7 7 €7 / A \ 7 €7 7 — / 7 €7 

2 V 



where we have represented in thick black (resp. red) line the initial (resp. mutated) boundary 
edge. 

In the above equations, the transformations b — > x (resp. v — > y) are precisely the two 
types of forward mutations of the A 2 T-system cluster algebra, obtained by applying the 
first (resp. second) line of ( 13. 1 1) . We may now turn to the proof of Theorem 13.41 

Proof. The formula is proved by induction under mutation. We start from the basic stair- 
case solution (13. 7J) , which may be put in the form ( 13 . 13j) . upon substituting A = A\A 2 and 
B = BiB 2 , and noting that k\ — 1, jx — j + k — 1, while ko = 1 and jo = j — k + 1. 

Starting from the expression (13.131) for the basic staircase boundary, we may apply iter- 
atively either of (13. 14j) or (13. 1 5|) to get to any other boundary (up to global translation), 
by simply substituting products of pairs of triangle matrices into the expression (I3.13p . 

We must however pay special attention to the extremal cases, namely when the mutation 
acts on the edge just before the upper extremity as in Figj6] (a), or just after the lower 
extremity as in Figj6] (b), of the projection of (1, j, k). 

The prefactor Ti,^,^ in (13. 1 3[) corresponds indeed to the bottom vertex of the upper 
extremity of the projection of (l,j,k) onto the boundary. Assuming as in Figj6] (a) that 
the edge just before the upper extremity of the projection of (l,j,k) is of e type, with 
value b at the bottom vertex as in eq. (13.141) . the mutation sends it to an /-type edge 
with bottom vertex value x, which becomes the new upper extremity of the projection of 
(l,j,k), replacing c. Noting that -E>i(x, c, v)j,\ = <5j,if, we see that the last multiplication 
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u 



u 




a A 



c 




c 



(1,j,k) 



X 



a 



x 



(1,j,k) 



(a) 



(b) 



FIGURE 6. The effect of a mutation next to the top (a) and bottom (b) of the pro- 
jection of (1, j, k). The new mutated edge is represented in red. The corresponding 
new lower vertex of the projection extremity is changed accordingly: c — > x (a) 
and a — > x (b). The extremal matrices Bi(x,c, u) (a) and Ai(a,x,u) (b) must be 
dropped from the expression for M(i,j), as the corresponding triangles lie outside 
of the new projection of (1, j, k). 

by Bi(x, c, u) amounts to replacing the global prefactor c by x, which is the desired change 
of T l jljfcl . Analogously, when the mutation acts on an edge of type e next to the bottom 
extremity of the projection as in Figj6](b), with bottom vertex value b as in (I3.14p . it sends 
it to an edge of type / with bottom value x, which becomes the new lower extremity of the 
projection of (1, j, k), replacing a. As Ai(a, x,u)ij = 5j t i, we may drop the contribution of 
this first triangle, and we recover ( I3.13p . This completes the proof of the Theorem. □ 

The case a = 2 needs no extra work, due to the following symmetry: 

Lemma 3.6. For any fixed boundary, with initial data of the form {T a j^ a ^\ a =\,2-,j&, we 
have 



Corollary 3.7. The solution of the A2 T- system for arbitrary boundary is for all a = 1,2 
a Laurent polynomial of the initial data with non-negative integer coefficients. 



4.1. T-system. We now consider the general A r T-system ( |l.mi.2p for j,k e Z, and 
a e I r . As before, this splits again into two independent systems for T a j^ with fixed value 
of a + j ' + k modulo 2, which we fix to be 0, without loss of generality. 

The indices {(a,j,k)} for a + j + k = modulo 2 run over r consecutive horizontal 
layers a = 1,2, ...,r of the centered cubic lattice, each of which is a square lattice, the 
vertices of the next layer lying at the vertical of the centers of the faces of the previous one. 
For technical reasons, we will also represent the extra bottom and top layers a = and 
a = r + 1, within which all values of T are fixed to 1, by the A r boundary condition ( II. 2p . 




Proof. The transformation a->r + l- aisa symmetry of ( ll.lHl.2p . 



□ 



4. The A r case 
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FIGURE 7. The basic staircase boundary for the A r T-system, here for r = 4, on 
which all vertex values T a jfi and T a j-\ t i are specified (blue dots). We have added 
an extra bottom and top layer where, according to the A r boundary, all vertex 
values are set to 1 (red dots). 

4.2. Boundaries. As before we start with the basic boundary {T a j t o, ^a,i-i,i}ae/ T . jez f° r 
a + j even. We may describe this boundary in 3D space as a succession of broken lines at 
constant j (represented in thick solid lines in FigJTJ) of the form: 

(4.1) ij = {(a,j,e aJ )} r a=1 , e aJ = a + jmod2 

We denote by A,B the vertical stacks of tetrahedrons depicted in FigJTJ respectively between 
(-2i-i and £ 2 i and £ 2 i and £ 2 i+i for all j £ Z. Each tetrahedron has two opposite (thick) 
edges of the form (a,j, e) — (a + 1 — e) and (a,j + 1, 1 — e) — (a + l,j + l,e). 

4.3. Solution for the basic staircase boundary. In Ref. [B], the system was solved for 
the basic staircase boundary in two steps. First one eliminates T a j^ for all a > 1 as: 

(4.2) T a jk = det (Txj-a+b fe+a+6 _ a _i) , a G I r , j,k G Z 

l<a,6<« 

which allows to concentrate on Tij^. 

The solution Txj^ wa s then expressed in terms of paths on a rooted target graph with 
time- dependent weights involving only the boundary values. More precisely, T lj)fc was 
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found to be equal to Tij+kfi times the partition function for paths starting from the root 
at time j — k and ending at the root at time j + k. It is best expressed in terms of the 
2r + 2 x 2r + 2 transfer matrix, encoding the step weights y = yi, y 2 , ■ 2/2r+i : 
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Defining the time-dependent weights 

,t— 1,1 -^a,t+l,0 



yi(t) 



2/2a+l(*) 



V2a(t) 



T, 



a+l,t 



i T a . 



l,t+l,0 



[a 



1,2, 



l,t+l,0 J- a+l,t,Q 1 a,t,l 1 a,t,0 J- a,t+l,l 

and the transfer matrix for steps from time t to t + 1: 

7(t,t + l)=7( yi (t),y 2 (t),...,y 2r+1 (t)) 

we have: 

Theorem 4.1 ([6]). The solution of the A r T -system $1.1W1.2\) for the basic staircase 
boundary reads for a = 1 : 

j+k-l 



T hjj k = T hj+kj0 ( n 



t=j-k 



4.4. Reduced transfer matrix. As before we note that the two-step transfer matrix 
T(y, y') = T(y)T(y / ) is again decomposable into two linear operators acting on two comple- 
mentary spaces of dimensions r+1, corresponding respectively to components (1, 3, 6, 7, 10, 11, 
and (2, 4, 5, 8, 9, ...). Explicitly, the operator acting on the first set of components reads: 



( y[ 


1 
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2/2 + 2/3 + 2/4 
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2/6 + 2/7 + 2/8 
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The corresponding reduced two-step transfer matrix is CP(t, t + 2) = CP(y(£),y(t + 1)). 
Theorem 14.11 turns into: 

k-i 

(4.3) T 1Jtk = T lj+M ( Yl 9{j - k + 2i,j -k + 2i + 2) S j 

4.5. Gauge transformation and rhombus/triangle decomposition. We may apply 
to CP any gauge transformation of the form 7{t,t + 2) = Lftit^t + 2)L^ t _ 1 2 and such that 
(L t )ij = (£t)j,i = <5j,i without altering the result (14. 3p . 
Here we choose: 

T — A- (~\ T ^2,i,1^3,t,0 rp rp ^2a,t,l72a+l,t,0 \ 

L>t — Cliagl^l, 1 x >t) Q, — — , 13^0, • • • , J-2a-l,tfi, — ™ > ■ ■ ■) 

•*3,i-l,l 2a+l,t-l,l 

One advantage of this choice is that CP(t, t + 2) only depends on values of T a j t k at times 
j = t, t + 1, t + 2. It is also justified a posteriori by the decomposition formulas below. 

We will now decompose the reduced two-step transfer matrix CP into a product of ele- 
mentary matrices, defined as follows. 

Definition 4.2. We define the following 2x2 elementary step matrices: 

(4.4) H(a,b,x)=(l 7(x,a,6)=U f 



Note that these generalize the horizontal and vertical step matrices of Definition 12. 2[ 
used in the Ai case. 

Definition 4.3. For any 2x2 matrix X and a G {1,2, ...,r} define the r + 1 x r + 1 







f 1«-1 


a -lx2 


0q— 1 xr- 


(4.5) 




02XQ-1 




02xi — a 






^ Or-axo-1 


Or— ax2 


If— a 



where l m denotes the m x m identity matrix and mxp the m x p matrix with zero entries. 

This gives rise to matrices H a ^ a+ i(a, b, u) and V a)a+ i(a, 2, y) using for X f)4.5f) the matrices 
(I4.4p . As an example, the triangle matrices A\, A 2 , B\, B 2 introduced in the A2 case (I3.10p 
may be identified with: 

Ai{a,b,u) = H lj2 (a,b,u) A 2 (b,u,v) = V 2 ^(b,u,v) 

Bt(a,b,u) = Vi )2 (l,a,6) B 2 (b,u,v) = H 2j3 (u,v,l) 

We finally define: 

Definition 4.4. We define the following r + 1 x r + 1 matrices: 

H at0l+ i(x,y,u) if a odd 
V a>a+ i(a, x, y) if a even 

V a)a+ i(a, x, y) if a odd 
H aj a+i(x,y,u) if a even 



(4.6) A a (a,x,y,u) 

(4.7) B a (a,x,y,u) 
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A A, A 2 A 3 ... B B, B 2 B 3 



FIGURE 8. The triangle decomposition of the stacks A,B of tetrahedrons on the 
basic staircase boundary. The triangles are colored either white or gray, and are 
grouped into pairs of distinct colors sharing an edge (rhombi). The H (resp. V) 
rhombi have the gray triangle pointing up (resp. down), and correspond to prod- 
ucts with the matrices H (resp. V) in the indicated manner. The three parameters 
are the values of T a j t k at the vertices of the gray triangle. 

and for sequences a = a , a%, a r+ \ and b = bo, b\, b r+ i, we define: 

A(a;b) = Ai(6 , a x , b x , a 2 )A 2 (b 1 , a 2 , b 2 , a 3 )...A r (b r ^x, a r , b r , a r+1 ) 
(4.8) 5(a;b) = B 1 {bQ,ai,b 1 ,a 2 )B 2 {b 1 ,a 2 ,b 2 ,a 3 )...B r {b r ^ 1 ,a r ,b r ,a r+1 ) 

Theorem 4.5. The following decomposition holds: 

?(t, t + 2)= A(a(t), a(t + 1)) B(a(t + 1), a(t + 2)) 

where a a (t) = T a j fia t , with e a ,t as in ( 14. ip . and in particular ao(t) = a r+ i(t) = 1, due to 
the A r boundary condition. 

Proof. By direct calculation. □ 

In analogy with the frise result, we may now interpret pictorially the decomposition 
result above. We interpret the matrices A(a(j — l),a(j)) and £>(a(j),a(j + 1)) as corre- 
sponding respectively to the stacks A,B of tetrahedrons in FigJTl Each tetrahedron may be 
decomposed in two ways into a pair of triangles (the thin solid lines in FigJT] correspond to 
one particular choice for each tetrahedron). 

Let us now pick the canonical choice indicated in FigJS] for both types A, B of stacks, 
namely that all squares are divided along their second diagonal. Then each term Aj or Bi 
in the products ( 14 .8p may be associated to a pair of triangles sharing an horizontal edge 
(which we call a rhombus, by a slight abuse). 

As illustrated in FigJHl we color the triangles in white or gray in such a way that rhombi 
are made of one triangle of each color. This gives rise to two types of rhombi H (resp. V), 
depending on whether the gray triangle is on top (resp. bottom), that are interchanged 
when going from stacks of type A to type B. In the rhombi, gray triangles play a particular 
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role: their vertex values carry the 3 parameters of respectively the matrices H or V in the 
product definition of A, or Bi (I4.6ti4.7p . Taking the product over the rhombi from bottom 
to top yields the formulas (14.81) . 

Each stack contains r — 1 tetrahedrons, hence there are 2 r_1 possible triangle/rhombus 
decompositions of each stack of type A (resp. B), each made of r rhombi. By their definition 
(14 . 5ti4 . 6tl477j) . the operators Ai and Aj (resp. Bi and Bj) commute as soon as \i — j\ > 1. 
So a given rhombus decomposition carries the information of whether Ai multiplies Aj_|_i 
from the left (like in FigJS]) or from the right (for the other choice of diagonal in the i-th 
tetrahedron from the bottom), and of what the three arguments of the H or V factors are 
(via the boundary values at the vertices of the gray triangle). 

The 2 r ~ l a priori distinct matrix products corresponding to these rhombus decomposi- 
tions turn out to be identical. This is a consequence of the following local commutation 
relations (for odd i): 

A i ^ 1 {x,a,b,a!)A i (b,a! ',b' \y) = Ai(a,a ,b' ,y)A i _ 1 (x,a,b,b') 
Bi-i(x,a,b,a')Bi(b,a',b',y) = Bi(a,a' ,b' ,y)B i ^ 1 (x,a,b,b') 

easily derived from the following: 

Lemma 4.6. 




Proof. By direct computation. □ 



Theorem 14.11 now turns into: 
(4.11) 

fc-i 

Tij,k = Tij +kfi ^ Y\_ ^( a — k + 2i), a(j — k + 2i + l))£?(a(j — k + 2i + l), a(j — k + 2i + 2))j 

i=0 
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FIGURE 9. A forward mutation of the boundary consists of a vertex value update 
y — >■ y 1 , with yy' = xz + au. It takes the 4 triangles in the back (black lines) of the 
octahedron (a) to the 4 triangles in the front (red lines). We have also represented 
(b) the corresponding Ai-Bi matrix identity, equivalent to a rhombus exchange, in 
the case of odd i. 



With the above-mentioned freedom of picking any of the 2 r_1 A and B product decompo- 
sitions for each transfer matrix A(a(t), a(t + 1)) and B(a.(t + 1), a(t + 2)) in (14. lip , we get 
2 fc ( r_1 ) equivalent expressions for Tij^. 

4.6. Mutations as rhombus exchange. Forward mutations of the boundary of the A r 
T-system may take place whenever five neighboring vertices form the back of an octahedron, 
say in positions k — 1), — 1, k), + 1, k), (i — k) and (i + 1, j, k), in which 
case the mutation fi it j replaces the back vertex (i, j, k — 1) with the front one (i, j, k + 1). It 
also induces the update of the boundary value y of the back vertex with y' = hu+xz for the 



front one (see FigJH] (a), with x = Tij_ hk , y = T idM _^ z = T ? 



T, 



u 



and y' = Tj Jifc+1 ). Backward mutations just correspond to the inverse process with y' — > y. 

Whenever a mutation is possible, using the freedom to pick rhombus decompositions, 
we will see in next section that we may always bring H and V rhombi in contact. The 
mutation corresponds then to interchanging the two rhombi: it is a forward mutation if V 
passes from the left to the right of H (see Fig. [9]), backward otherwise. 

To check that the mutation is faithfully represented by the matrices, we use the following: 

Lemma 4.7. 




(4.12) 



Proof. By direct calculation. 



V(b, x, y) H(y, z, u) = H(x, y', u) V(b, y\ z), yy' = xz+bu 



□ 



The forward mutation //jj is implemented by simply switching the two rhombi in the 
decomposition of the corresponding adjacent stacks. In the case i odd, this corresponds to 
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the identity depicted in Figj9](b): 

Ai(b, x, y, u) Bi(b, y, z, u) = B^b, x, y', u) A^b, y', z, u) yy' = xz + bu, 

a direct consequence of Lemma 14.71 When i is even, the mutation /ij j corresponds to the 
same identity read in the opposite direction (but still corresponds to a V passing from the 
left to the right of an H). 

We conclude that the forward mutations of the boundary may be iteratively imple- 
mented on the formula (14.111) by simply (i) picking the relevant rhombus decomposition 
among all equivalent ones namely bring a V to the left of an H with the vertex (i,j,kij) 
of the boundary in the center (ii) switching the corresponding matrices and Bi in the 
product transfer matrix. This will be made very precise in the next section. 

4.7. Arbitrary boundaries and the 6 Vertex model. In this section, we give a bisec- 
tion between boundary strips of the A r T-system and configurations of the 6- Vertex (6V) 
model on infinite strips of width r — 1. In the 6V picture, mutations correspond to the 
reversal of all spins around square faces whose spins form an oriented loop. 

4.7.1. Local configurations and the six vertices. Recall that the 6V model is a statistical 
model defined on the square lattice, whose configurations consist of an orientation (spin) 
on each edge, in such a way that the following "ice rule" is satisfied: at each vertex of the 
lattice there are exactly two entering and two outgoing edges. 

The idea of the bijection is very simple. As we noted before, the boundary strips may 
be decomposed into broken lines £j at time j, j G Z, as in (14.11) . each made of a succession 
of r — 1 edges of type e or /. As discovered in the case of A 2 , there are 6 possible 
transitions between an edge of type e or / of the form e = (a, j, k) — (a + k + 1) or 
/ = (a, j, k) — (a + l,j, k — 1), at time j, connecting vertices in the layers a and a + 1, 
and an edge of either type at time j + 1, also connecting vertices in the layers a and a + 1. 
The two corresponding edges give rise to either tetrahedrons A,B if they are of opposite 
types, and parallelograms C,D,E,F otherwise. Moreover, the 4 vertices included in the two 
edges also define two horizontal edges respectively within the layers a and a + 1. These 
may also be of only two types, say e' and /', of the form: e' = (a, j, k) — (a,j + 1, k + 1) 
and /' = (a, j, k) - (a, j + l,k- 1). 

These 6 local configurations A,B,C,D,E may be represented as squares in projection 
onto a (j, a) plane, with vertical edges labeled e or / and horizontal edges labeled e' or 
/'. In this projection, a complete boundary is simply a strip of squares, infinite in the j 
direction, and of finite width in the a direction, with some compatible edge assignments. 
The corresponding configurations of the 6V model are on the dual of this strip, namely 
with a vertex in the middle of each square. 

To the 6 configurations of each square, we associate bijectively the 6 vertices of the 6V 
model as indicated in FigJTOj namely we pick the horizontal edge orientation to be to the 
right (resp. left) for an e type (resp. / type) edge, while the vertical edge is oriented 
upward (resp. downward) for an e' type (resp. /' type) edge. 
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e' f f e' f e' 




f e' f e' f e' 



A B C D E F 

FIGURE 10. The projections in (a, j) plane of the six local configurations of the 
boundary of the A r T-system and the corresponding local arrow configurations of 
the dual 6V model. We have indicated for each configuration the corresponding 
tetrahedron (A,B) or parallelogram (C,D,E,F). 




FIGURE 11. The 6V configuration corresponding to the basic staircase configu- 
ration of FigfT] for r = 4. 

Another direct way of connecting the initial boundary vertices to the 6V configuration 
is to view the latter as coding the value of k a ,j as the label of the face (j, a) in the 6V 
configuration as follows. Note first that values of k differ by ±1 in neighboring faces. We 
impose the following "Ampere rule" that the value of k on the right of an arrow is smaller 
than that on the left. This fixes all the values of k a j up to a global translation. 

4.7.2. The basic staircase boundary as a 6 V configuration. We now consider the basic stair- 
case boundary introduced in Section 14.21 It is an infinite horizontal strip of square lattice 
with height r — 1, with an alternance of e, f, e, f, ... vertical edges and e', f, e', f ... 
horizontal edges, namely a checkerboard of configurations A and B of tetrahedrons, with 
vertices such that k a j G {0, 1}. The dual 6V configuration is simply made of arrows 
that alternate in all directions (it is a fully antiferromagnetic groundstate configuration). 
In particular the bottom and top horizontal rows of vertical edges alternate along the j 
direction. 

4.7.3. Mutations as loop reversals. As illustrated in Figj9j a mutation of the boundary at a 
given vertex amounts to transforming the 4 edges sharing this vertex, by simply permuting 
the two vertical ones (i.e. belonging to the plane j =const.) and the two horizontal ones 
(i.e. belonging to the plane a =const.). This is also true if a = 1 or a = r, with the 
condition that the bottom (resp. top) vertex which belongs to the layer a = (resp. 




FIGURE 12. A sequence of mutations applied on a length 2 portion of the bound- 
ary strip for A2 in the 6V formulation. The green (resp. red) dots mark the places 
where a forward (resp. backward) mutation can be applied. The spin of the top 
and bottom edge sequences are conserved, equal to —1 and 1 here. 

a = r + 1) has an attached value 1. This situation was already encountered in the case 
r = 1 (see Lemma 12. 6ft and r = 2 (see Lemma [3. 51) . 

Such a mutation can only take place if the permuted edges are of respective following 
types, clockwise around the vertex from the top: e, e', /, /' for a forward mutation, and 
/, /', e, e' for a backward mutation, and the two are transformed into one-another. 

Once reformulated via the 6V bijection, and for 1 < a < r, this simply means that a 
forward (resp. backward) mutation can take place only around faces whose adjacent edges 
form a clockwise (resp. counterclockwise) oriented loop, and the mutation simply reverses 
the direction of all 4 arrows. When a = 1 (resp. a = r), the "loop" is open, i.e. is only 
formed of 3 edges with the bottom (resp. top) edge of the loop missing. In the case of A\, 
we are only left with open loops with both the top and bottom horizontal edges missing. In 
both cases, the corresponding mutation reverses the orientations of the 2 or 3 edges around 
the loop. By a slight abuse of language we still call "loops" these open loops, and still refer 
to the edge reversal as "loop reversal" . 

Starting from the basic configuration of FigJTTj we may therefore generate all others by 
successive elementary loop reversals. For illustration, we display in FigJT2]a sequence of mu- 
tations applied to a length 2 portion of boundary in the A2 case. The central configuration 
corresponds to the basic staircase boundary. 

All our 6V boundaries are generated by iterated loop reversals on the 6V basic staircase. 
However, the basic staircase has the particular property that its top and bottom boundary 
vertical edges alternate between up and down. Loop reversals involving these will create 
zones of successive up or down edges, with the property that the total spin (i.e. the total 
number of up minus down edges) remains 0. We deduce the following: 

Lemma 4.8. The boundaries of the A r T-system correspond to the configurations of the 6V 
model on an infinite strip of width r — 1, such that there exist two times j m i n and j ma x £ ^ 
with the same parity such that: 

• (i) The portions of 6V boundary for j < j min and for j > j max are identical to those 
corresponding to the basic staircase. 
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FIGURE 13. The correspondence between 6V configurations and the 
square/triangle decomposition of the dual lattice. The diagonals that split squares 
into two triangles are the diagonal reflection symmetry axes of the vertex configu- 
rations (dashed lines). Vertices without diagonal reflection symmetry correspond 
to squares. The gray (resp. white) triangles or squares correspond to horizon- 
tal arrows pointing toward (resp. from) the central vertex, and vertical arrows 
pointing from (resp. toward) the central vertex. 

• (ii) The total spin for the portion for j min < j < j max of 6V boundary is zero for 
both top and bottom vertical edges. 

Finally, the 6V boundaries must also carry the information of the initial data: there is 
one initial value per vertex of the original boundary, hence these may be represented as face 
labels for the 6V configuration (including the bottom and top faces with only 3 bounding 
edges) . 

4.8. From general boundaries to networks. We now wish to associate transfer matrices 
to the 6V configurations described above. 

4.8.1. 6V and square/triangle decomposition. To define the transfer matrices associated to 
6V configurations, we first associate to each 6V configuration a square/triangle decompo- 
sition of the dual square lattice via the correspondence of FigfT3j Note that each edge of 
the square lattice is adjacent to two polygons (triangle or square) of opposite colors. 

4.8.2. Slice transfer matrix. Any given 6V boundary is decomposed into vertical slices of 
unit width, made of a vertical sequence of r — 1 vertices. To each such slice we associate 
anr + lxr + 1 matrix T. The matrix T is defined recursively from the bottom to the top 
of the slice according to its 6V configurations as follows. 

For illustration and help, we have represented in FigJTH a typical slice configuration (a) 
in square/triangle decomposition and (b) in 6V form. We first consider the square/triangle 
decomposition of the 6V slice at hand, the vertices of which carry the initial data, say 
ax, a,2, a r from bottom to top on the left and b\, b%, b r on the right. We further split 
all the squares into pairs of triangles by drawing their second diagonal. Note that horizontal 
edges are adjacent to two triangles of opposite colors. By a slight abuse of language, we call 
again rhombi such pairs of triangles. With this definition, we are left with two unpaired 
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(a) (b) (c) (d) 

FIGURE 14. A typical vertical slice in square/triangle decomposition (a) and the 
associated 6V configuration (b) in the Aj case. The squares in (a) are further 
cut into two triangles by drawing the second diagonal. The configurations (a)-(b) 
determine the H, V structure of the transfer matrix (c) where green (resp. red) 
segments between lines + l correspond to -f^j+i (resp. V^+i) matrices. The 
relative position of segments in consecutive horizontal slices corresponds to the 
order of multiplication of the associated matrices. We also indicate the network 
(d) corresponding to the complete matrix product. The labels dj, b\ are now face 
labels of the network. 

triangles with edges a r — b r and a\ — b\ respectively. This is repaired by adding the missing 
triangle of the opposite color on both top and bottom of the slice, with a third vertex 
carrying the value a r+ \ = b r+ \ = 1 and = bo = 1. As they all carry the same value 1, it 
will be convenient to identify all the vertices of the bottom layer as well as all those of the 
top one. 

We label the lines supporting the horizontal edges of the 6V and the two extra rows of 
vertices by 1,2, ...,r + 1, from bottom to top. We denote by Tj the matrix corresponding 
to the partial slice extending from lines 1 to i. 

We start by defining 7\ = Hi^(cli, &i, x) if the bottom triangle is white and Tj = 
^1,2(1,01,61) if it is gray, where x is the top vertex value of the rhombus containing the 
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bottom triangle. Let us consider the vertices i and i + 1, for i = 1,2, ...,r — 1. If the 6V 
vertical edge between line i and i + 1 points up (resp. down), then T i+1 is obtained by 
multiplying Tj by an operator V iji+ i (resp. H ii+ i) of Eqs. (I4.9ti4.10j) . Moreover the order 
of multiplication is from the left (resp. right) if the diagonal of the square is the first, 
connecting the SW and NE corners (resp. second, connecting NW to SE corners). This 
gives: 





(4.13) 



1+1 1,1+1 i 



Hu+Ji 



TM 



,i+i 



TV, 



j+i 



TV, 



TH 



i+1 



The three arguments of the V i>i+ i (resp. H i>i+1 ) operators are (a^-i, a iy bi) (resp. (a i; b i: u i+ i)), 
where (resp. Ui+i) is the third vertex label of the gray triangle with edge a* — bi. This 
defines T = T r uniquely. 



Example 4.9. Let us consider the case of Fig[T4\ (a)-(c). 
obtained recursively as follows: 



The total matrix T = T 7 is 



T 1 = V 1 , 2 {l,a 1 ,b 1 ) T 2 
T5 = H$fi (05 , 65 , &g) T4 Tq 



V 2 ,3(ai,a 2 ,b 2 )T 1 T 3 
T 5 V 6 j(bQ,a 7 ,b 7 ) T 7 



T 2 V 3A (b 2 ,a 3 ,b 3 ) 
H7,8(a 7 ,b7,l)T 6 



T4 — T 3 ^4,5(04, 64, 65 



4.8.3. Networks. As pointed out earlier, the operators with respective indices i, i + 1 and 
j, j + 1 commute as soon as j > i + 1 or j + 1 < i. A useful way of representing the matrices 
Vi j + i and H i i+ i is as networks. We start with the representation: 



i+1 




(4.14) 



H iii+1 (a,b,u) 



V iji+1 (x,a,b) 



where the left vertices i, i + 1 are connected to the right vertices i, i + 1 via weighted edges 
i — > j corresponding to the nonzero entries (i, j) of the matrices H, V of f)4.4p . More 
precisely, a diagonal entry 1 with indices (j, j) is coded by a thick black horizontal edge 
connecting the left vertex j to the right vertex j. The other nontrivial entries {j,k) are 
coded by thick green (resp. red) edges for H (resp. V), connecting the left vertex j to 
the right vertex k. The actual values of these nontrivial entries are coded by the labels of 
the faces of the network x, a, b, u in (I4.14p . Implicitly, the full representation of Hi^ + \ and 
Vi i+i asr + lxr + 1 matrices also involves thick horizontal black edges connecting the left 



THE SOLUTION OF THE A r T-SYSTEM FOR ARBITRARY BOUNDARY 



31 



vertices 1,2, — l,z + 2,...,r + lto their right counterparts, corresponding to the block 
form of (14. 5p . and which we omitted in (I4.14p for simplicity. 

The multiplication of two matrices is naturally coded by the concatenation of their 
networks, by identifying the right vertices of the network of the first matrix to the left 
ones of that of the second. The result is in general a rectangle of height r + 1, and width 
possibly smaller than the number of matrices multiplied, as some horizontal black lines may 
be freely removed to gain space. We denote by N(M) the network associated to the matrix 
M. For illustration, we have represented in FigJTH (d) the network N(T) corresponding to 
the matrix T of Example 14.91 As another example, the "mutation" Lemma 14.71 reads in 
network language: 

u u 

zs = 

b b 

The networks provide us with natural weighted path models. Indeed, the matrix entry 
T m p is interpreted as the partition function for paths from left to right on the network of T, 
starting from the left vertex m and ending at the right vertex p, weighted by the product 
of weights of the edges visited, namely: 

(4.15) T m , p = II w ^ 

paths m— >p edges e 
on JV(T) visited 

4.9. General solution as path model on networks. We now give the general expression 
for the solution T a j ^ of the A r T-system for arbitrary boundary conditions. 

4.9.1. The case of As shown above, any boundary condition is coded by a config- 
uration of the 6V model on an infinite strip of width r — 1 with the properties (i)-(ii) of 
Lemma 14.81 

Definition 4.10. The projection of (l,j,k) on the boundary is the portion of boundary 
between the broken lines £j and £j l , containing the bottom vertices (1, Jo, ko) o,nd (1, jx, ki), 
repectively such that jo — ko = j — k and j\ + k\ = j + k with jo maximal and j\ minimal. 

To this projection we naturally associate the corresponding truncated configuration C of 
the 6V model on a rectangle of height r — 1 between the planes j = jo and j = ji, and with 
face labels given by the original vertex values of the boundary. Let Tj 0J1 (C) denote the 
product of vertical slice transfer matrices from the slice (jo, jo + 1) to the slice (ji — l,j%). 
Then we have: 

Theorem 4.11. The solution of the A r T-system with arbitrary fixed boundary reads for 
a = 1: 




(4.16) 

in the above notations. 



Tij,k — ^iji,fei {fjoji (^))i,i 
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Proof. This is proved by induction under mutation. First, (14. 16ft is satisfied for C = Co, 
corresponding to the truncated basic staircase boundary. Indeed, we may rewrite (14. lip as 

| J B(a{j - k + 1 + 2i), a(j - k + 2 + 2i))A(a(j - k + 2 + 2i), a(j - k + 3 + 2i)) 



v i=0 



1.1 



by use of A (a, = Sj^ and -B(a, a')^^ = Sj^a^/ai. We see that this boils down to (14.161) 
with k = ki = 1, ji = j + k — 1, and j = j — k + 1. As explained above, any mutation 
H may be implemented by an elementary loop reversal on the 6V configurations, and only 
mutations within the projection of k) affect the value of Txj^- 

Let us assume ( 14.161) holds for some 6V configuration with face labels C. We wish 
to apply a mutation /i, i.e. form the configuration //(C), identical to C except for the 
reversed elementary loop and the updated face label. Assume this is a forward mutation, 
i.e. the corresponding loop is oriented clockwise in C, and assume for definiteness that it 
occurs between vertical slices s and s + 1, with respective slice transfer matrices T^ s ' and 
T( s+1 \ and between horizontal lines i and i + 1. Due to the rules (14.131) . examining the 
triangle decompositions of the four adjacent vertices to the loop, we find the following four 
possibilities: 















\ 




Recalling that the choice of diagonal in the white squares is arbitrary, we may bring all 
situations to the first one, by use of the identity (I4.10p within the relevant slice transfer 
matrices. Applying the mutation to this case amounts to the transformation of Figj9] (b). 
We now have the following four possibilities for the environment of the reversed loop: 



















) 












To recover the correct slice transfer matrices, we must flip the diagonals of the gray squares, 
by applying (14. 91) . Taking the product of transfer matrices over all slices, we finally get the 
"mutated" matrix Tj 0J1 (/i(C)). We may repeat the same argument for backward mutations, 
while exchanging the roles of white and gray triangles, and for mutations at vertices with 
a = 1 and a = r with the obvious changes. 

As before, we must however consider separately the case when the mutation occurs in 
the lower left or lower right square face of the truncated 6V configuration. In these cases 
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indeed, the projection of (l,j,k) is modified and we must drop (resp. insert) a first or 
last slice transfer matrix when the mutation is forward (resp. backward). As before, 
dropping/inserting a last slice transfer matrix implements the correct change of prefactor 
Txj u kx that makes it agree with the new projection. In all cases, (14.161) follows for //(C). 
This completes the proof of the Theorem. □ 

4.9.2. Non-intersecting paths on networks and the case ofT a jj,. Let Nj 0} j 1 (C) denote the 
network associated to the transfer matrix T J0J1 (C). Then T\ j t k/Ti j 1 ^ is interpreted as the 
partition function for weighted paths on the network N^^C), that start at the left vertex 
1 and end at the right vertex 1, namely: 



(4.17) e n 

1 iji,fci 



w(e) 



paths p: 1— )-l edges e 
on j (C) visited 



where w(e) stands for the weight of the edge e. 

Note that in the case of the basic staircase Co, this path interpretation is different from 
that of Ref. [6]. Nevertheless, we may like in Ref. [6] interpret the determinant identity 
(14. 2p that relates T a j^ to Tij^ in terms of non-intersecting paths, now on networks. 

Let us pick an arbitrary boundary C, and rewrite the determinant formula ( 14. 21) as: 

Tg,j,k _ ( -?lJ-a+b,A;+a+6-a-l 



Ul<b<a T l,Mb)M{b) l<a.6<« V T 1Jl{b)M{b ) 

where we denote by (1, ki(b)) the common rightmost bottom vertex of the projections 
onto the boundary C of the vertices (1, j '• — a + b, k + a + b — a — 1), a = 1, 2, a. Similarly 
let us denote by (1, jo(a), fco(a)) the common leftmost bottom vertex of the projections onto 
the boundary C of the vertices (1, j — a + b, k + a + b — a — 1), b — 1, 2, a. 

We denote by (j, x) the vertex of iV(C) that lies between the slices j — 1, j and j,j + 1 
at height x, namely the common exit and entry point x of respectively the network for the 
slice j — 1, j and that for j,j + 1. For instance the left vertex x of N a ^ b (G) has coordinates 
(a,x), while the right vertex of same height has coordinates (b, x). 

We finally have: 

Theorem 4.12. The solution T a j^ of the A r T -system for arbitrary boundary is equal 
t° rii<6<a ^di(b)M(b) times the partition function for families of a non-intersecting paths 
on the network iV(C), starting at the vertices (jo(a),l), a = 1,2,..., a and ending at the 
vertices (ji(b), 1), b = 1, 2, a. 

Proof. In view of (I4.16P and its reformulation in terms of paths on networks ( I4.17p . we may 
interpret the (a, b) term in the determinant (I4.18P as the partition function for paths from 
the bottom left to the bottom right vertex on the network Nj ^j 1 ^(Q). The Theorem 
follows from the application of the Lindstrom-Gessel-Viennot theorem [16] [12]. □ 

Corollary 4.13. The solution T a j^ of the A r T -system for arbitrary boundary condi- 
tion T a j )ka . = a a j (a e I r ,j G Z) is a positive Laurent polynomial of its initial data 

( a a,j)a£l r ;jeZ- 
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(a) (b) 

FIGURE 15. A view in perspective of a particular boundary (a) for the ^3 T- 
system (blue vertices and thick solid black edges). We have indicated the three 
horizontal layers a = 1, 2, 3, and the j, k coordinates at a = 1. We have represented 
the vertex P = (2, 3, 3) and the four points of the a = 1 layer Q\ = (1, 3, 2), Q% = 
(1,4, 3), Q3 = (1,2, 3), Q4 = (1,3,4) involved in the determinant formula (|4.19p . 
The projections of the latter onto the boundary are X = (1,0,1), Y = Q\ = 
(1, 3, 2) and Z = (1, 6, 1). We also represent the same information (b) in projection 
onto the (k,j) plane, with layers a = 1,2,3 represented respectively in black, red, 
blue. We have also indicated in dashed lines the mutations leading to this boundary 
from the basic staircase one. 



4.10. An example. Let us consider the case r = 3 and the boundary depicted in FigJT5| 
with the A 3 boundary condition that T j^ = T 4ji fc = 1 for all j, k. We wish to compute 
T 2j3i 3 in terms of the indicated boundary values (blue vertices): 



21,0,1 


= ai 


21,1,0 


= &i 


21,2,1 


= Ci 


21,3,2 


= di 


21,4,1 


= ei 


21,5,0 


= /l 


21,6,1 


= 91 


21,0,0 


= a 2 


21,1,1 


= b 2 


21,2,2 


= c 2 


21,3,1 


= d 2 


21,4,0 


= e 2 


21,5,1 


= h 


21,6,2 


= 92 


2~3,0,1 


= as 


23,1,0 


= h 


21,2,1 


= c 3 


2~3,3,0 


= d 3 


21,4,1 


= e 3 


2~3,5,0 


= h 


23,6,1 


= 93 



We have the determinant formula 



21,3,2 


21,4,3 




2qi 


2q 2 


21,2,3 


21,3,4 




2q 3 


2q 4 



which involves the 4 points Qi,Q2,Q3,Qa of the a = 1 plane. In turn, the projections of 
the Qi onto the boundary read: X = (l,j (2), fc„(2)) = (1,0,1), Y = (1, j (l), k (l)) = 
(1,^(1)^1(1)) = (1,3,2) and Z = (l,ji(2), h(2)) = (1,6,1). 

Next, we construct the 6V configuration with face labels of Fig [16] (a) corresponding to 
the boundary. The simplest way is to start from the basic staircase and apply to it the 
3 mutations that bring it to the boundary under consideration (see the dashed lines in 
Fig{15] (b)). We also construct the triangle decomposition of Fig[T6] (b), by drawing the 
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FIGURE 16. The 6V version of the boundary of Figll5l(a) with its face labels and 
the corresponding triangle decomposition of the dual (b). The associated network 
(c) is represented with an indication of the successive slices (vertical dashed lines) 
and the corresponding face labels. The truncation (d) of the network is the only 
information needed to get the result. The partition function for paths from X to 
Z that do not visit Y is a weighted sum over the four paths (e),(f),(g),(h). 



diagonal reflection symmetry axes of the vertices and a second diagonal in the remaining 
squares. This allows to build the network of Fig{16] (c), by following the above procedure. 
In particular the edge-weights of the network are given by the surrounding face labels 
according to the rules f!4.14p . 

The Lindstrom-Gessel-Viennot theorem expresses T 2t3 ^/ (T lt3j2 T ltGtl ) as the partition func- 
tion for pairs of non-intersecting paths on this network that start at X, Y and end at Y, Z. 
This is simply the partition function for a single path from X to Z that does not visit Y. 
To compute the latter, we may restrict ourselves to the truncation of the network depicted 
in Fig{16] (d), where the edges keep their original weights from (c). The partition function 
is therefore a weighted sum over the four paths (e),(f),(g),(h), and the final formula follows: 

(4 20) T 2,3,3 = _J_ + c 3 d 3 + c 2 eie 3 + c 2 f 2 

71,3,2^1,6,1 d 3 gi d 2 d 3 g x d x d 2 e 2 g x d x e 2 g x 



5. Application: Q-system from the T-system solutions 

In this section, we compare the restriction of our T-system solution to the case of the 
A T Q-system whose solution was worked out in [I]. We find alternative path models for 
describing the general solutions. 
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5.1. Q system. The Q-system for A r is obtained from the T-system by "forgetting" about 
the variable j. This can be done in various ways, here we adopt the following: we impose 
that the solution and the initial data be periodic, namely that T , Q! ,i+2 ) fc — ^aj.fc f° r a U a , J, fo- 
lia this case, the quantities R a ^ = T aja+ k[2],k are easily seen to satisfy the A r Q-system: 



In [I], the boundary data for (15. ip were shown to be in bijection with Motzkin paths of 
length r — 1, m = (mi,m 2 , ...,m r ), and to take the form x m = {R a ,m a i Ra,m a +i}a=v The 
explicit solution R a ^ for each such data was expressed as follows. First, for each Motzkin 
path m one constructs an explicit rooted oriented graph r m with 2r + 2 vertices, and with 
weighted edges encoded in a 2r + 2 x 2r + 2 transfer matrix T m . The root corresponds 
to the row and column index 1 of the transfer matrix. The matrix T m is constructed in 
such a way that any oriented edge % — > j of T m going away from the root ( "ascent" ) has 
a trivial weight 1, while any oriented edge % — > j pointing toward the root ("descent") has 
some-non-trivial weight t x ^ (m), the latter being a Laurent monomial of the initial data. 
The main result of JI] for a = 1 reads: 

(5.2) R\^ n+ini = Ri tmi ((/ — T m ) ^ 

where the notation X\pi denotes the coefficient of t n in X. In other words, the quantity 
Ri. n+mt/ Ri,mi is the partition function for paths from and to the root on T m , with exactly 
n descents. 

Example 5.1. For the case of A 3 , the transfer matrix T 2j i 5 o for the Motzkin path (2, 1,0) 
reads [1]: 



(5.1) 



Ra,k+lRa,k-l 
Ro,k 



R(x,k + Ra+l,kRa-l,k. 
Rr+l,k = 1 



(5.3) 



/o 

tyi 








1 

1 

ty 2 1 

ty 3 

£3/3,1 tyi o 



*3/4,l ^2/4,2 








1 




ty& 
o 



o o 

o o 





1 1 





ty 7 









1 



0/ 



where: 



Vi = 



Rl,3 



,2/2 = 




,2/3 = 




,2/4 = 




,2/5 = 




-^2,1-^1,2-^1,3 



-^3,0-^2, 1-^2,2-^1,3 
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The general solution for R a , n was shown in [I] to be expressible as partition function for 
families of a strongly non-intersecting weighted paths on T m , via a generalization of the 
Lindstrom- Gessel- Viennot Theorem. 

5.2. Network formulation from T-system. As mentioned above, the Q-system restric- 
tion amounts to having a periodic solution of the T-system, with period 2 in the variable 
j. For such a solution, boundaries are much simpler, as they are also 2-periodic in the 
j direction. However, constructing such a boundary by iterated mutations on the basic 
staircase (which has the desired periodicity) involves repeating each mutation an infinite 
number of times (with periodicity 2 in the j direction). So stricto sensu these boundaries 
are not covered by our previous solution. But for each fixed value of (a,j, k) only a finite 
portion of the boundary is needed to express T a So we may apply a finite number of 
mutations for each case and get the correct result, and just formally complete the boundary 
by periodicity. We may therefore apply the results of Section H] here. 

In the periodic case, the boundary values must read T a ^ )ka0 and T a ^y. al for j = 0, 1 
mod 2 respectively. Note that k a ^ — k a ^ = ±1. Defining m a = Mm(k at0 , k a ^), m a + 1 = 
Max(fc Q , j o, k a ,i), we find that m is the relevant Motzkin path for describing the Q-system 
boundary data. 

According to our results, each periodic boundary may be viewed as a configuration of 
the 6V model but now with a periodicity of 2 in the direction of the strip. In other words, 
we have a configuration of the 6V model on a cylinder of perimeter 2 and height r — 1. 
Moreover, as we started from a configuration with alternating vertical arrows (spin with 
two edges) on both the bottom and the top of the strip, any mutated configuration has one 
vertical edge pointing up and one pointing down on the upper and lower boundaries of the 
cylinder. We have the following: 

Lemma 5.2. The 6V configurations on a cylinder of perimeter 2 and height r — 1 with 
alternating edge orientations on top and bottom are in bisection with Motzkin paths of length 
r — 1, with mi = or 1. 

Proof. The bijection goes as follows. We start from a Motzkin path m = (m^mj, ...,m r ), 
with mi = or 1. Then we have three possible situations for each of the r — 1 steps of the 
Motzkin path: Xi = m^+i — = 0, —1, or 1, and two possible values for y$ = % + rrii = 0, 1 
modulo 2, for i — 1, 2, r — 1. We have the following dictionary between the six possible 
pairs (xi,yi) and the six vertices of the 6V model. 




(0,1) (0,0) (1,1) (1,0) (-1,1) (-1,0) 
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2 




2,1 



3,0 



1,2 



1 



(a) 



(b) 



FIGURE 17. The periodic 6V configuration (a) associated to the Motzkin path 
(2,1,0) for r = 3. The corresponding boundary vertices (a,j,k) may be read off 
the configuration: the third index k is determined by the Ampere rule and by 
m\ = 2, and is indicated in each face. We have also sketched the associated rhom- 
bus/triangle decomposition (b) of the two corresponding slices, with the assigned 
boundary values. 

More simply, Xj is the "algebraic sum" of the two horizontal edges (0 if they are opposed, 
1 if they both point to the right, —1 if they both point to the left) and yi is determined 
only by the bottom vertical edge (1 if it points down, if it points up). 

Let us arrange the vertices associated to {(xj, ?/j)}i=i on a single column, from bottom 
to top. Note that the rules above make the orientations of edges compatible, so we can 
identify the bottom vertical edge of the vertex i + 1 with the top vertical edge of the vertex 
i. There is a unique way to complete this configuration into one on the 2 x r — 1 cylinder. 
We must indeed add a second column of vertices which are identical to the previous ones, 
up to reversal of all vertical arrows (the unique solution respecting the ice rule), in order to 
satisfy both the horizontal periodic boundary condition and the alternating one on top and 
bottom. Note that mi = or 1 determines whether the bottom left vertical edge points 
down or up. □ 

Example 5.3. We consider the case r = 3 and the Motzkin path (2, 1,0). We have X{ = 
= 1 for i = 1,2. The 6V configuration is represented in Fig[T7\ (a). We have 
indicated the values of k a j in the faces of the configuration, corresponding to a = 1,2,3 
and j = 0,1. These are not to be mistaken for the usual face labels of the 6V configuration, 
which are the assigned boundary data a a j for the corresponding vertices: T a j^ a = a a j . 

To each Motzkin path m = (mi, mi, ui r ), we may now associate a configuration of 
the 6V model on a cylinder of perimeter 2 and height r — 1 via Lemma 15.21 by using 
the shifted Motzkin path m' = m — 2(p,p, ...,p), where p = p-^]. The corresponding 
values of k a j are uniquely determined by the Ampere rule and by m a = Min(A; a)0 , k a! %), 
m a + 1 = Max(& a) o, fc Q ,i)- The solutions for of Theorem 14. 1 1 1 and for T a of Theorem 
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14.121 involve only two types of slice transfer matrices, due to the periodicity, namely that 
for all slices of the form [2a, 2a + 1] and that for [2a — 1, 2a]. These are coded by the two 
columns of the 6V configuration, which have the same horizontal edge orientations and 
opposite vertical ones. Using the construction of Sect. 14.8. 2\ the slice transfer matrices 
may be constructed inductively as follows, directly from the Motzkin path m: 

We start from mi. If it is even, the first slice transfer matrix has T\ = Hi^, if it is odd 
it has T\ = Vi t 2- Then, having constructed T i} say with a last factor H i)i+ \ (resp. V^+i), 
we have three possibilities (according to the value of x^ = m i+1 — mi): 

• (i) Xi = 0: then T i+1 = V iyi+1 Ti = TjV^+i (resp. T i+X = H iA+1 Ti = T*if M+1 ). 

• (ii) Xi = l: then T i+X = TiH iA+1 (resp. T i+1 = V^+iTj). 

• (iii) Xi = -1: then T i+1 = H iyi+1 Ti (resp. T i+1 = TjV^+i). 

The arguments of the matrices are, as usual, the boundary data at the vertices of the 
gray triangles, with the vertex indices of the form (a,j, k a j) where k a j are determined by 
m. This gives a matrix U m = T r for each Motzkin path m. The second slice is treated 
analogously. Due to the fact that the its 6V configuration is identical to the first up to 
reversal of all vertical arrow, it is easy to write the corresponding slice transfer matrix 
U m in terms of U m . Indeed, the rhombus/triangle decomposition of the second slice is 
the reflection of the first w.r.t. a vertical axis, and with all colors of triangles inverted. 
We therefore have to interchange H -h- V, and to reverse the order of the factors. More 
precisely, let * be the involutive anti-automorphism ((ab)* = b*a*, (a*)* = a) such that 
V* i+1 = H i>i+ i, then we have U m = (U m )*. 

Example 5.4. In the case of the Motzkin path (2, 1, 0) forr = 3, we easily read the transfer 
matrices on the rhombus/triangle decomposition of Fig[F3\ (b): 



(5.5) 



2,1,0 



U- 



2,1,0 



-^3,4(^3,0, -^3,1, 1)^2,3(^2,1) -^2,2, -^3,1, 

\ 






Hl t 2\Rl,2i #1,3, R'. 



2,2 J 



/ 1 

R 2 ,2 




Rl,2 





R 2 ,l 


Rl,3 

Ra,i 


Rl,3 

Ri, 2 R3,i 


Rl,3 

\rTJ 


R13R2.2 

Rl,2 


R2,2 
R2,l 


Rl,iR2,2 


^2,2^?3,1 



^3,0 1 

Rs,i/ 



^1,2(1, Rl,3, #1,2)^2,3 (Rl,2, R2,2, #2,1)^3,4(^2,1, #3,1, R. 



3,0j 



/ Ri,3 

1 Rl,2 




V 



R 2 ;. 



Ra, 



Rl 1 2R2.1 
R 2 ,2 

R 2 ,i 






R2,lR3,0 

Ri, 2 R3,i 

R2A,R3,0 

Rsa 
Rs,o 




R3,0 » 
fll,2 

Rs,o 
R 2 ,i 

Rs.o 
1 / 
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Let us apply the result of Theorem 14.111 with ko = ni\ + 1 = k\ and k = n + mi, 
jo — j — n + 1, ji — j ' + n — 1. Defining e = 1 — {k mod 2) and 9 = j\ mod 2, we get: 



Rl,n 



+ nii 



Ti t e,k — 



(U m U, 



\n-l 



\n— 1 



1.1 



1.1 

R\, mi (U m U n 



1.1 



where we have used (U m ) liX = 5 liX and (U m ) Xy i = ^.l^i.i-e.&i-i/^i.e^i = 5x,i- R i,m 1 /-Ri, mi +i- 
Comparing this with (I5.2p . we deduce the 

Theorem 5.5. Lei m be a Motzkin path of length r — 1, T m the 2r + 2 x 2r + 2 transfer 
matrix of the Q -system solution of Ref. [1], and £7 m , U m as above. Then we have an 
identity between "resolvents" : 



((I-T m )- 1 ) ll = Ul-tUnU, 



i.i 



Note that the matrices U m and U m have size r + 1 x r + 1. We may however view the 
product U m U m as the transfer matrix of a weighted graph G m with 2r + 2 vertices labeled 
1, 2, (r + 1), 1', 2', (r + 1)', defined as follows. We interpret the matrix element (U m )ij 
as coding the weights of the oriented edge i — > j 1 of G m while (U m )i,j codes the weight 
of the oriented edge i' — > j. Alternatively, we may form the 2x2 block transfer matrix 

Ur, 



e 



m 

u m 



j for G m and note that — \/t0 m ) x j 



{I-tU m U n 



i.i 



Noting that the non-zero elements of U m have the same indices as those of the transpose 
of U m , we see that G m has doubly oriented edges only, with specific weights for each 
orientation. From the network construction, all these weights are Laurent monomials of 
the initial data. 

Example 5.6. For the Motzkin path m = (2,1,0) of r = 3, we may use the matrices 
(15. 5p . Without altering the resolvent, we may gauge-transform the matrices U m and U m 
with invertible matrices R,L with Lj t i = L 1; j = 5jy. V m = L~ l U m R and V m = i? _1 f/ m L. 
For the choices: 



t j- /i R1.2R2A -^2,1-^3.0 

L = diag 1, — , — . 



R 



2,2 



R 



R 



we find that 



Vr, 



( 1 

V2 
2/3,1 





2/3 
2/4 



3,1 




2/5 



3,0) R — diag(l, i? lj2 , i?2,i, -^3,1/ 



0\ 






Vrr 



\2/4,l 2/4,2 2/6 1/ 



/2/1 1 1 1\ 

111 

11 

\0 y 7 J 



in terms of the parameters of (15. 4p . One checks directly the statement of Theorem \5 . 51 with 
the matrix (15. 3p and by computing the rational fraction ( (I — tV m V m ) 



1.1 
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Finally, turning to R a , n for a > 1, we have a simpler picture than in jl], as the quantity 
-R«,n+mi/(-Ra,mi) a is expressed directly as the partition function for a non-intersecting paths 
on the associated network, without having to generalize the Lindstrom-Gessel-Viennot 
Theorem. 

6. Conclusion 

In this paper we have presented an explicit solution of the A r T-system in terms of 
arbitrary boundary data. This solution is expressed in terms of partition functions of 
weighted paths on some particular networks, determined by the boundary. 

We have briefly described the connection of the T-system to a particular cluster algebra. 
In particular, the sets of boundary data we have considered here form only a subset of the 
clusters in this cluster algebra. What happens is that the form of the cluster mutations 
can change in general from the equation fll.ip . in some sense, the equation itself evolves, 
leading to clusters of a different kind. Nevertheless, the positivity conjecture seems to hold 
for these other clusters as well. It would be extremely interesting to probe whether these 
other clusters also have a description in terms of networks, that would make the Laurent 
positivity property manifest, like in the cases studied in this paper. 

Another possible direction of generalization concerns non-commutative cluster algebras. 
In [TJ, the cluster algebra for the non-commutative A\ Q-system was introduced, and its 
positivity proved by use of a non-commutative weighted path model. We have checked that 
it can be reformulated as a non- commutative network model in the spirit of the present 
paper, but with transfer matrices with entries in a non-commutative algebra. We hope to 
report on this direction in a later publication. 
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